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Abstract 

A nonlocal gravity model is considered which does not assume the existence of a 
new dimensional parameter in the action and includes a function /(D -1 /?), with □ the 
d'Alembertian operator. Using a reconstruction procedure for the local scalar-tensor 
formulation of this model, a function / is obtained for which the model exhibits power- 
law solutions with the Hubble parameter H = n/t, for any value of the constant n. 
For generic n — namely except for a few special values which are characterized and also 
specifically studied — the corresponding function / is a sum of exponential functions. 
Corresponding power-law solutions are found explicitly. Also the case is solved in all 
detail of a function / such that the model contains both de Sitter and power-law solu- 
tions. 



1 Introduction 

A large number of observational data [U [2j [3j HJ [5j [6] supports that the expansion of the 
Universe is accelerating now and has been doing so for a considerable part of its existence. 
The assumption that General Relativity (GR) is the correct theory of gravity leads to the 
conclusion that about seventy percents of the energy density of the present day Universe 
should be a smoothly distributed, slowly varying cosmic fluid called dark energy |7J. If, 
alternatively, one wants to describe the accelerating expansion without this fluid, the gravity 
equation must be modified (for recent reviews, see, e.g., [HEUdO]). There has been a lot of 
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discussion on which side of the Einstein equations should be modified, either by adding the 
cosmological constant, or anything similar, to the matter side of the equation, or by modifying 
the gravity side of the equation [U]. Note that some of these modified gravity models, for 
example, F(R) gravity models, can be transformed into general relativity with additional 
scalar fields by a suitable conformal transformation of the metric [TOT [12]. 

Higher-derivative corrections to the Einstein-Hilbert action are being actively studied in 
the context of quantum gravity (as one of the first papers we can mention [13], see also [UJ and 
references therein). A nonlocal gravity theory obtained by taking into account quantum effects 
has been proposed in [15]. In a way, this development goes back to string/M-theory, which 
is usually considered as a possible theory to encompass all known fundamental interactions, 
including gravity. The appearance of nonlocality within string field theory provides therefore 
a good motivation for studying nonlocal cosmological models. 

The majority of available nonlocal cosmological models explicitly include a function of 
the d'Alembertian operator, □, and either define a nonlocal modified gravity [TJ)J [TBI El EIB1 
HH1 [201 ED 1221 1231 EH 123 [23 I2U or contain a nonlocal scalar field, minimally coupled to 
gravity [28] . 

In this paper, we consider a nonlocal gravity model that includes a function of the d -1 
operator and does not assume the existence of a new dimensional parameter in the action [15] . 
For this model, a technique for choosing the distortion function /(d -1 ), so as to fit an arbitrary 
expansion history has been derived in [20]. The considered nonlocal model has a local scalar- 
tensor formulation [T7] . 

It has been shown in [17] that a theory of this kind, being consistent with Solar Sys- 
tem tests, may actually encompass all known sequence of the Universe history: inflation, 
radiation/matter dominance and a subsequent and final dark epoch. In [T9] the ensuing cos- 
mologies corresponding to the four basic epochs: radiation dominance, matter dominance, 
acceleration, and a general scaling law have been studied for nonlocal models involving, in 
particular, an exponential form for the function f(ip). An explicit mechanism to screen the 
cosmological constant in nonlocal gravity was discussed in [22 [ [23 | I2T]. 

The cosmological model often includes some function which cannot be deduced from the 
fundamental theory, usually one can only assume some property of this function, for example, 
that the function should be a polynomial or more generally an analytic function. Then the 
question arises, why some specific form of the function /(D -1 ) is chosen and not a different 
one, what is the physical motivation for the final choice. To answer these questions one has to 
resort to the reconstruction procedure itself and to understand in depth the expansion history 
of the Universe in order to match both with care. In the framework of this local formulation, 
a reconstruction procedure has been proposed [19J in terms of functions of the scale factor a. 
The reconstruction method considered in this paper was first proposed in [25] with the aim 
to obtain models with de Sitter solutions and it has been generalized in [27] ■ This method 
uses convenient functions of the cosmic time t. 

In the present paper, we develop in some detail the reconstruction procedure and use it to 
find models with power-law solutions. We obtain functions / which correspond to power-law 
solutions, e.g., solutions for which the Hubble parameter is given by H = n/t, where n is a 
real number. Also, we will find functions / for which the model has both the de Sitter and 
power-law solutions. 
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It should be noted that the reconstruction procedure is widely used in cosmology, in 
particular, in models with minimally coupled scalar fields [29J [30j [3TJ [32] , also nonminimally 
coupled scalar [33] or Yang-Mills [34] fields, F(R) and Gauss-Bonnet gravity models [32l [35] . 
and F(T) models, with T being the torsion scalar [36] . 

The paper is organized as follows. In Sect. 2, we consider the action for a general class 
of nonlocal gravity models and derive the corresponding equations of motion in the Jordan 
frame. In Sect. 3 we propose the algorithm which allows to find f{ip) and get the model 
with the given cosmological evolution. In Sect. 4, we prove that the model with f(ip) being 
a sum of exponential functions can indeed lead to power-law solutions. This is extended in 
Sect. 5 where we find sufficient conditions on f(ip), under which the model has power-law 
solutions, and present these solutions in the explicit form. In Sect. 6 a function f(ip) is 
considered which corresponds to a model with both de Sitter and power-law solutions, what 
is a quite interesting, somehow unexpected finding. In Sect. 7, we consider a few special cases 
of power-law solutions and get for each case the corresponding function f(ip). The last section 
is devoted to conclusion. 



2 Scalar-tensor gravity model 

We consider a class of nonlocal gravities, whose action is given by 

S = J d A xV=g [R (1 + fia-'R)) - 2A] + £ mattcr | , (1) 

where k 2 = 8ttG = 8tt/M P i 2 , the Planck mass being M Pi = G~ 1/2 = 1.2 x 10 19 GeV, while 
/ is a different iable function that characterizes the nature of nonlocality, with D _1 being the 
inverse of the d'Alembertian operator, A is the cosmological constant, and £ ma tter the matter 
Lagrangian. We use the signature (— , +, +,+),<? being the determinant of the metric tensor 
g^y. Recall the covariant d'Alembertian for a scalar field, which reads 

where V M is the covariant derivative. 

Introducing two scalar fields, ij) and £, one can rewrite action (CD) as follows 

S = J d A x^g~ I JL + /(V>) + - - 2A] + £ matter j . (2) 

By varying action (j2J) over £, we get Oip = R. Substituting ip = D _1 i? into action (J2|), one 
reobtains action ([!]). Action (J2J) may thus be regarded as a local action. 

By varying this action with respect to £ and ip, respectively, one obtains the field equations 

□V> = R, (3) 

□£ = f'ty)R, (4) 
where the prime denotes derivative with respect to ip. 
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Varying action (J2J) with respect to the metric tensor g^ v yields 
l -g, v [m + d p ^ - 2(A + □*)] -R^V - \ {d^d v ijj + d^d^+V^ = -K 2 T mflv , (5) 
where \1> = 1 + f(if>) + £, the energy-momentum tensor of matter T mpu is defined as 

rp _ 2 (5 (\/— g£ m attcr) /„x 

The system of equations (jSJ) does not include the function ■0, but only f(ip) and f'(ip), together 
with time derivatives of ip. Note that f(ip) can be determined up to a constant, because one 
can add a constant to f(ip) and subtract the same constant to £, without changing any of the 
equations. 

To consider cosmological solutions we assume a spatially flat Friedmann-Lemaitre-Ro- 
bertson-Walker (FLRW) universe, with the interval, 

ds 2 = -dt 2 + a 2 (t)5 ij dx i dx i , (7) 

and consider scalar fields which depend on time only. In the FLRW metric, the system of 
Eqs. (J3J) — (JSD reduces to 

3H 2 m = -^-3^ + A + K 2 p m , (8) 
(2H + 3H 2 ^j V = -H>-V- 2H4? + A - K 2 P m , (9) 
| = -3Hi-Q{H + 2H 2 )f^) 1 (10) 
# = - 3HtP - 6 (H + 2H 2 ) , (11) 



where differentiation with respect to time t is denoted by a dot, and the Hubble parameter is 
H = a/ a. We assume that matter is a perfect fluid with T m00 = p m and T m y = P m gij- The 
equation of state (EoS) reads 

p m = -3H(P m + p m ). (12) 
Adding up (jSJ) and we obtain the following second order linear differential equation for 

* + 5if# + (2H + 6# 2 ) * - 2A + /t 2 (P m - p m ) = 0. (13) 

In the case when A = and matter is absent, Eq. ffl3|) has a solution \I/(t) =0. This solution 
is not interesting, because, as one can see from Eq. (jSj), the corresponding ij) and £ are both 
constants. 



3 Nonlocal gravitational models with a given cosmolog- 
ical evolution 

Let us now consider the question, how to get the function f(ip) that corresponds to a solution 
with a given Hubble parameter H(t)l If the EoS parameter w m is known, then for any given 
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H(t) one can integrate Eq. (|T2|) and get p m (t), after which one can solve Eq. (fT3|l and obtain 
^(t). The knowledge of H(t) also allows to integrate ([3]) and get ip{t). 
Substituting 

i{t) = *(t) - /(vo - 1 

into Eq. ffTUl) and using the function t(ifj), we get a linear differential equation for f(ip), namely 

/' (^ip 2 - 12 (H + 2H 2 ^j f{$) = ^ + 3F^ . (14) 

It can be more suitable for the calculation to rewrite Eq. (fT4l) as follows: 

sip- 12(H + 2H 2 )s = $ + 3H4f, (15) 

where s(t) = f'(ip(t)). From this first order differential equation one can obtain the function 
s(t), and, after that, get f'(ip) and f(ip). 

Note that Eq. (fT4"j) is a necessary condition for the model to have the solutions in the 
given form. Substituting the obtained function f(ip) into Eq. ([8]) or Eq. fl9]), one can check 
the existence of solutions of the given form. Note that this method not only allows to construct 
the corresponding f(ip), but it also gives a suitable way to obtain £(t) and ip(t) for a given 
H{t). In the following Sections we will illustrate this result for some physically important 
examples of possible H(t) behaviors. 



4 Power-law solutions 

In [25] we have found the general form of the function f(ip) for which the model can have de 
Sitter solutions. Here, we explore power-law solutions, in other words, we assume that the 
Hubble parameter is of the form H = n/t. 

We consider matter with the EoS parameter w m = P m / p m being a constant, not equal to 
— 1. For H = n/t, Eq. (1121) has the following general solution: 

p m (t) = p t~ 3n ^ +1 \ (16) 

where po is an arbitrary constant. Eq. ({TBI) has the following general solution: 

• For n 7^ — 1 and n ^ -1/3 

A nnK 2 (w — 1 \t 2 -3(i+w m )n 

m/ = c x r 2n + c 2 t l ~ 3n + £ -t 2 - , PoK (17) 

(n + l)(3n + 1) (3nw m — l)(n + 3nw m — 2) 



For n 



^ = Cl t 2 + C 2 t 4 - At 2 ( ln(t) + -)- Pok 2 (w 1) 5+3Wm (lg) 



2/ 3(l + w m )(3w m + l 
For n = —1/3, 

* 2 = Cl t 2 + + ht 2 L( t ) -D- ^f^;' 1 !,/ ^- (19) 

2 V 4 / {w m + l)(3w m + 7) 
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where C\ and Ci are arbitrary constants. w m is chosen so that (3nw m — l)(n + 3nu> m — 2) be 
not equal to zero. 

Inserting H = n/t into Eq. (ITT]) and assuming n / 1/3 and n 7^ 1/2, we find the solution 
for ip{t) to be 

m = ^. 6 J^-Jl ia (l), (20) 



3n — 1 \to 

where and to are arbitrary constants. In all formulae we assume that t > 0, so, to > 0. 
Also, we get that the general solution of ffiTl) is 

= C 3 r 1/2 + C 4 , at n=i (21) 

V>(*) = ^(ln(t)) 2 + C ( 3ln(t) + C ( 4 , at n = ± (22) 

where C3, C4, C3, and C4 are arbitrary constants. We assume that ^(i) is not a constant, so 

To get /(VO in terms of elementary functions for any n, except for the above-mentioned 
particular values, it is useful to set ipi — in ( 12"U]) and to express the time variable t in terms 
of ip, as 

(l-3n) , 

t = toe 6 ^ 2 "- 1 )^. (23) 
We obtain the following general solution to (|14p : 

U2 _ 2j.2-3ra— 3mu m 
/q^) = c (l-3w)V»/(3n(2n-l)) _ FO' 1 ^0 e (3n(l+«; m )-2)(3n-l)i/)/(6n(2n-l)) + 



6n(l + n) 3n(n — 2 + 3nw m ) 

(W - l)t Q 2 "Cl (3n-l)V>/(3(2»-l)) , D (3n-l)V/(3n(2n-l)) , n 

2(2n-l) 



(24) 



where .Di and -D2 are integration constants. Note that C\ is an arbitrary constant as well. 
To get this result we assume that n 7^ 1/2, n 7^ 1/3, n 7^ —1/3, and 

Let us consider the case n = — 1. Substituting \I/(t), given by (|T8|) . into Eq. (fl4"]) . we get 
the following general solution: 

/"( a 2j9(5+3ui m ) 

hW = ^fa 4 ** + D ie ^ - A (is ln(to) +3 + 4^) - e^^^ + D 2 . 

3 54 w m + 1 

(25) 

Note that this form of f{ip) is similar to a function for which there exist de Sitter solutions at 
w m = 1/3 (see Sect. 6). At the same time, is a sum of exponential functions provided 

A = 0. 

For n = —1/3, we use (ITS]) and get 

9 1 n ^2j.3+U! m 

/ 2 ty,) = t C X fl*e**lS _ £ t 2 Ae 6^/5 _ gPO^ e | (3 +^ + ^^12^/5 + ^ 

5 4 3w m + 7 

Note that ^(V*) coincides with /o(V') for n = —1/3, so can use the formula (jUj) for n = —1/3 
as well. 
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The cases n = 1/3 and n = 1/2, where the function has a special form, different from 
( 120]) . will be considered in Sect. 7. 

We see that for all values of n, except when n = —1, n = 1/3, and n = 1/2, f(ip) is equal 
to either an exponential function or to a sum of exponential functions plus a constant. This is 
the general form of the function f(ip), which can have power-law solutions in the model with 
the cosmological constant and an ideal perfect fluid. We have not proven that this model has 
power-law solutions, what we really proved is that other models do not have such solutions 
(with ipi = 0). 



5 Power- law solutions for a given f(ip) 

The full system of equations is fl8|)-f lT2|) . We solve equations (1I0]) - (JI2"]) and the sum of equa- 
tions (jSJ) and (jUJ). Thus, to get the solution it is enough to substitute the functions obtained 
into Eq. flS}. 

After this substitution, we get that terms which are proportional to A and po cancel and 
that Eq. (JHJ) takes the following form: 

6tl (3n ~ 1) {2n-l)nD 1 t- 6n = 0, D x = 0. (26) 

Note that the case n = 1/2 has been excluded (in this case the function if>(t) has another 
form). Thus, we get that the model with 

fty) = f\e a ^ + + f 3 e a ^ + / 4 , (27) 

where fi and 014 are constants, has solutions with H = n/t, where n is an arbitrary real 
number, but n 7^ —1, n 7^ 1/3 and n 7^ 1/2. As we have noted in Sect. 2, the function f{ip) 
is defined up to a constant, so we can put J4 = without loss of generality. 

The method here considered allows not only to get the suitable function f(ip), but also 
to obtain solutions in explicit form. Indeed, we have found that the model considered, with 
f(ip) given by ( |27|) . has the solution 



H(t) 



n 

T 

6n(2n-l), ft 



At2 / 1 

" + ~, T, r , Tl 7^ 

(n + l)(37i + l)' 3 
6 + -At 2 ( ln(t) - 1 



(2f 



n = — 



where 



2 V 4 / ' 3 ' 

^(t) = *(t)-/(^)-l, 

2n , n +1 -3n Po^ m - l)^-3(l + ^)n 



e = dr zn + c 2 t 



(3nw m — l)(n + 3nw m — 2) 
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the constants being subject to the following conditions: 

~M 2 1 - 3n 



6n(l + n) 3n(2n-l)' 

p ^r 3n ~ 3 "" m (3n(l + ™ m )-2)(3n-l) 



/a = - o ; " , o — 7, 



3n{n - 2 + 3ra/; m ) ' 6n(2n - 1) ' (29) 

~ = (w - l)tp 2n Ci 3n-l 
/s 2(2n-l) ' " 3 3(2n-l)' 

h = o. 

Thus, for the model with the given f{ip) and A, the values of the parameters C\ and to are fixed 
and there is a connection between a\ and a 3 . The values of a 2 and f 2 define conditions on w m , 
which is not equal to —1, and Pq. Note that the solutions include an arbitrary parameter C 2 - 
Recall also that n ^ — 1, n / 1/3 and n ^ 1/2, and we demand that n — 2 + 3ray m 7^ and 
w m 7^—1- As a consequence of these results, it follows that the model without matter can 
develop power-law solutions only if f 2 = 0, thus f(ip) is either an exponential function or a 
sum of two exponentials. 



6 Models with both power-law and de Sitter solutions 
6.1 The function f(ip) for the model with de Sitter solutions 



In [25] we obtained that the model admits the de Sitter solution H = H , where H is a 
nonzero constant, provided that 

f d sW = Cl e^ + c 2 e^ - f k e^+W , for w m ? - \ , (30) 

3(1 + 3w m )iio <J 

fM) = Cl e^ + c 2 e^ + e^, for w m = - I , (31) 



ml V v ) ' m 3 



and where c\ and c 2 are arbitrary constants, p~o is the initial value of the matter energy density. 
Note that we include in this model only matter with the same EoS parameter w m , but that 
the initial values of the matter energy density, po and po, can be actually different. 

In [25], a de Sitter solution in its most general form has been found for f(ip) being an 
exponential function. In particular, we got that equations (fTTj) - (}T3"j) . which do not depend on 
the function f(ip), have the following solutions: 



8 



Hit) = Ho, 

p m (t) = j5oe -3(i+^)^o(t-*o) ) 
ij;(t) = -4H (t-t ), 



9(t) 



T 



K 2 po(w m - 1) 

3Hj)W m (l + 3w 

ii 

K PO 3H (t-t ) ( 

Ho [ 

T + 4/C PO c _2ff (t-to) 



-3H (w m +l)(t-t ) 



w m 7^ 0, w m 7^ -1/3, 



(32) 



0. 



3# 



(t-t ), w m = -l/3, 



where 



T = c 3 e' 3Ho{t - ia) + 4 Cl e- 2Ho{t ~' to) + 



A 

3#r 



(33) 



C3 and to being arbitrary constants. 

Substituting into Eq. fl8]) the solution ([32]) and f (t) = — — 1, where f(ip) is 
specified by (l30l) or ( 13~ll . one gets 



l2c 2 Hle- m ° {t - lo) = 0, 



c 2 = 0. 



Note that we obtain the same results for all values of w m . Thus, the model has de Sitter 
solutions if 



K 2 p 



3(1 + 3w m )H 2 



D 3(u> m +l)V>/4 



, for w m ^ 



-^fl \^\e^\ for 



AH, 







1 

3 



(34) 
(35) 



with arbitrary c\. 

Let us check for the possibility that the model, with some given function f(i/j), has both de 
Sitter and power- law solutions. Such possibility exists only in the case w m 7^ — 1/3. Indeed, 
for w m = —1/3 de Sitter solutions exist only for the function 



M„ 



(36) 



with 02 = 1/2, whereas a power-law solution exists only when f3 2 = 4/9. Therefore, a model 
of this type cannot have both de Sitter and power-law solutions. We thus conclude that the 
function f(ijj) should be either a single exponential function or a sum of two exponential 
functions. Let us now check these two possibilities separately. 
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6.2 A single exponential f(ip) 
Let us consider the model with 

/(VO = /oe w , (37) 

where fo an are constants. This case where f(i/j) is an exponential function has been mostly 
investigated [TT1 HH1 HH1 E21 [231 [2H 1251 [261 EZ] - Particular de Sitter and power-law solutions 
were found in [17] and, after that, de Sitter solutions have been studied in [231 [23 [26] and 
power-law solutions in [23J [37], but always separately. Using mainly results in [251 EI]> our 
aim here will be to construct a model which simultaneously exhibits both kinds of solutions: 
de Sitter and power-law. 

If matter is absent, then the model can have de Sitter solutions only in the case when 

It is easy to see that 0:3 7^ 1/2, for any n, so that the model can have power-law solutions 
only if «i = 1/2. In this case, 

1 V57 , \ 

"="4*12- (38 > 
Thus, we come to the conclusion that for p(t) = the model has both de Sitter and power-law 
solutions if and only if /3 — 1/2 and A 7^ 0. Note that fo can be arbitrary and does not get 
fixed by the value of A. 

In the case, when matter is present we get the condition 

a 2 = -(w m + l), (39) 

that gives 

4 

n = — -. (40) 

3(3 -w m ) 1 ; 

Note that the equation oti = ct2 with the n obtained has as only solution w m = — 1. As a 
consequence, these models can admit power-law solutions if A = 0. It has been found in [37J 
that, in the case A = 0, power-law solutions exist only if either n = l/(3w; m ) or if n — 1 and 
u> m = —1/3. Comparing these condition with f T40l) . we obtain that the model has both de 
Sitter and power-law solutions if and only if n = l/(3w m ) and w m = 3/5. In this case, we get 
P = 6/5 and n = 5/9. 

Thus, we have reached the conclusion that the model with a single exponential function 
has both de Sitter and power-law solutions provided that either p(t) = or A = (but not 
both at once). In the first case, (3 = 1/2. In the second one, (3 = 6/5 and w m = 3/5. 

6.3 A sum of two exponential functions 

Let us consider de Sitter and power-law solutions for 

/(VO = he^ + he^- (41) 
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We have shown that de Sitter solutions for such f(ip) exist only when matter is present. From 
( 154"|) . we obtain 

ft = ~, fa = l(w m +l). (42) 

Note that w m ^ — 1/3, hence (3i ^ (3 2 . 

The de Sitter solution obtained is given by ( 1521 . where 

~ 3(l+3 Wm )ff 2 - 4 

Cl = /l, PO = 2 J2) = ^P2 - 1, 

to and C3 being arbitrary parameters and H Q an arbitrary nonzero constant. 

Let us consider possible power-law solutions in this model. Because a 3 7^ 1/2, for any n, 
we need set a\ = 1/2, thus the value of n is fixed by ( 1551) . From a 2 = f3 2 we then connect n 
and w m by (j40l . to get 



25±3V57 , 

'"■^TW (43) 

The sign "— " corresponds to w m ~ —0.517, what looks realistic for cosmological models, 
whereas the sign "+" corresponds to w m ~ 4.517. 

In summary, we have found that the model with f(i/j) given by (1411 has both de Sitter 
and power-law solutions (a new and quite interesting case) if and only if both A and p(t) are 
not equal to zero and w m is given by ( 1431 . 



7 Special cases of power-law solutions 

In previous sections, when looking for power-law solutions, we have ignored some special 
values of n. For completion of our search we will consider those values of n in the present 
Section. Being precise, formula (|27|) is not valid for some (we call them here special) values 
of n. We discuss in what follows such special values and look for the corresponding functions 
/('••)• 

We first consider solutions with H = l/(2t). In this case the function i[)(t) is given by 
( 12~T|) . Substituting this ij)(t) and the function \l> in the form (fTT|) . into (1151) . we obtain: 



s{t) 



Ci 



4(3w m - 1) 



-p K t 



2.1-3?« m /2 



16 



C 3 Vt 3(3w m -2)C 3 ru ~ " 15C 3 
where D\ is an arbitrary constant. Relation (12 ip can be rewritten as 

Using this relation, we integrate ( 14"4"j) and obtain 



At 5/2 + D u 



(44) 



Ci 

a 



2 ^ 



C4C1 

2 ^r - °\ 



4ACf 



+ 



4/t 2 po 



15 - C 4 ) 4 3(3w m - 2) VV> - C, 



1— 3w n 



■ (45) 



The function f(ip) is defined up to a constant, so one can set the integration constant D 2 = 0. 
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By direct substitution into (E}, the following condition on the constants results: 

C 2 = C 3J Di + 6^. (46) 

^3 

This relation fixes the arbitrary parameter in the expression (|17p for the function ^(t). Thus, 
we get the function f(ip) which includes four arbitrary parameters Ci, C3, C4 and D\. In 
particular, we have found that if f(ip) is an arbitrary quadratic polynomial, then the model 
has solutions with H = l/(2t) if and only if A = and p(t) = 0. 

Let us consider solutions with H = 1 / (3t) . Substituting the functions ip(t), given by (I22l) . 
and \l/(t) into Eq. (fl~5l) . we obtain the following solution: 

^= [21n(t)W (i^3-l + 21n(t)]Af + 

9p n 2 [(w m - l)(21n(t) + 3C 3 ) + 2\_ , Wm 2d [3 + 2 ln(t) + 3C 3 ] ~\ ' ' ' } 
+ 3^5 1 W + 

Note that for w m = 1 the term which is proportional to po can be actually considered as part 
of the integration constant D\. 

Equation ([H]) can be rewritten in the following way 

= -- ~ si?jij -3H^/ + A + K 2 p m , (48) 

therefore, the knowledge of the functions H(t), ^(t), s(t), ip(t), and p m (t) is enough in order 
to check the existence of solutions in the given form. One does not need to use the explicit 
form of f{ip). This property is useful in the case when the form of f(ip) is too complicated, 
for example, when it is given by quadratures. From ( 1481) we obtain the condition: D\ = 6C2. 
After that, we use the relation between t and if) 



t = e { - 3d3±J)/2 , where J = yj 9C% + I2ip - 12C 4 , 

to get f'{ip)- Note that the function ip(t) has a minimum, what means that we get two inverse 
functions. 

Integrating the function thus obtained over ip, we get 

= -^^el^- 1 ^ (e^"^ +Ei(l,±hw m -l)J 



3w m — 5 \ \ 2 



+ (e*i + Ei (l, ±£) \ + (e ±J + Ei (1, T J)) (49) 

+ ^^(3^1 + 4^-4^) +D 2 , 

where D2 is a constant and the function Ei(l, —x) = — Ei(x), where the exponential integral, 
for all non-zero values of x, is defined as follows: 

e z 

Ei (.r) = / — dz. 
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For positive values of x the integral has to be understood in terms of the Cauchy principal 
value, owing to the singularity of the integrand at zero. This example clearly shows that the 
method here considered allows to obtain the functions f{ip) even in the case when they are not 
elementary functions, what is a bonus. We also get the power-law solution in explicit form. 



8 Conclusion 

In this paper we have extended the reconstruction procedure for the scalar-tensor model 
proposed in [17J, which is a convenient local formulation of the popular nonlocal model [15] . 
We have explicitly shown that the models with f{ip) given by fl271) with (129]) contain power-law 
solutions with H = n/t. The solutions obtained include an arbitrary parameter, so that what 
we actually find is a full one-parameter family of solutions. This result has been obtained for 
an arbitrary n, except for the values n = —1, n = 1/3, and n = 1/2. For those values of n 
the function f{ip) has a different form, which has been studied as well. 

Another result of our work is that some models with f(if>) being a simple exponential 
function have both de Sitter and power-law solutions, but this occurs if and only if either 
the cosmological constant A or the matter energy density p m (t) are equal to zero. At the 
same time, the model with a single exponential f(ip), which does not include neither the 
cosmological constant nor matter contributions, can have either de Sitter solutions or power- 
law ones, but never both of them simultaneously. These results characterize in a clear way 
all different possibilities. On top of these findings we have also demonstrated that there 
are models with both de Sitter and power-law solutions which include both a cosmological 
constant and matter terms. The function f{ip) for this kind of models is namely the sum of 
two exponential functions, specified by ( l4TT) and ( 1421) . 
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